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Abstract 

This paper joins some concepts from Mechanics, Partial Differential 
Equations and Control Theory in order to solve bi-time optimization prob¬ 
lems related to stress tensor in plastic deformations. The main goal is to 
analyze some optimal control problems constrained by the equilibrium 
equations of the stress tensor in perfect plastic plane medium. As conse¬ 
quence of this approach, a natural split of the constraints arises, leading 
to integrability conditions and changes a classical variational problem into 
an optimal control one. The final outcomes confirm all the expectations 
related to the physical features of plastic deformations phenomenon. 
Keywords: multitime maximum principle, complete integrability condi¬ 
tions, non-linear PDE system, perfect plastic plane medium. 

Mathematics Subject Classification (2000): 70H06, 70H30, 70Q05, 

49J20. 

1 Introduction 

Over the last few years, a research team from University Politehnica of Bucharest, 
supervised by professor dr. C. Udriste, has proved some very interesting the¬ 
oretical facts related to multitime optimal control m-m. [la-iiB]). Lately, 
in order to overcome the theoretical dimension and to achieve some practical 
confirmation, the team focused on applying these results for meaningful prob¬ 
lems from different scientific areas: differential geometry, statistics, mechanics 
etc. m ,[I3],[I9]). This paper is the consequence of these efforts, proving the 
applicability and utility of a general multitime maximum principle and also 


I 


emphasizing remarkable particularities of the specific mechanical problem ap¬ 
proached here. 

Section 2 analyzes the complete integrability conditions for a quasi-linear 
(non-elementary) PDE system with two variables (a bi-time) and a multi¬ 
dimensional undetermined map (2-sheet state variable). The result obtained 
here is applied in order to describe the integrability context for the PDEs defin¬ 
ing the 2-dimensional stress tensor. The major outcome of the section consists 
in emphasizing a natural split of a non-linear PDE system, over state gradients 
and a natural insertion of control variables, leading to manageable integrability 
conditions. Section 3 describes optimization processes constrained by non-linear 
PDEs, a bi-time maximum principle. The main result is adapted for constrained 
variational problems. Moreover, this section points out the natural transforma¬ 
tion of variational processes to optimal control ones, via the canonical controls 
resulted from integrability requirements. Finally, Section 4 applies the theo¬ 
retical facts from the first two sections in order to describe the solution of a 
variational process involving the stress tensor in perfect plastic plane medium. 


2 Complete integrability conditions 
for plane quasi-linear PDE systems 

The main goal of this Section is to phrase the complete integrability conditions 
for the PDE system describing the stress tensor in perfect plastic medium. Basi¬ 
cally, this starts with finding a pertinent description for integrability conditions 
for an abstract, general problem. Three important consequences shall derive 
from our attempt: 1) there is a natural split of an arbitrary quasi-linear PDE 
system, separating state gradients one from each other, generating manageable 
integrability conditions, 2) this natural split is determined by a natural insertion 
of control variables and 3) any variational problem constrained by an arbitrary 
linear PDE system may be naturally rephrased as an optimal control problem 
via the control variables mentioned above. 


2.1 One state-variable plane PDE systems 

Let D be a compact subset of R^, with global coordinates t = and let 

E = 3D denote its boundary, li A = : V, x R x ^ ^)a,p=T 2 

B = {B°‘ : D X i? X ^ R)a=T2 denote some 2x2, respectively 2x1 tensors, 
they define the following one state quasi-linear PDE system my- 


(1) A{t,x{t),u{t))gTadx{t) = B{t,x{t),u{t)), 


where, for the sake of simplicity, t stands for and denotes the bi-time 

variable^ x represents the state variable (unknown differentiable two-sheet, in¬ 
volved in the PDEs via its gradient) and u{t) = {u^{t))i.^YN denotes the control 
variable (parameter for the PDE system). Written explicitly, the PDEs are: 


( 1 ') 


I A^^{t,x{t),u{t))^ -j- A^‘^{t,x(t),u{t))^ = B^{t,x{t),u{t))-, 

I f^T f^T 

y A^^{t,x(t),u{t))-^ P A^‘^{t,x{t),u{t))-^ = B‘^{t,x{t),u{t)). 
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Suppose X = is a solution for the foregoing PDE system, defined 

implicitly by ^{t,x{t)) = 0. 

Applying the Implicit Functions Theorem, it follows 


dx 


9$ 

dx 


a = 1,2. 


Substituting these in (!') leads to a linear PDE 


9$ 


,12 


A^^—+A^^—+B^— = 0 


dB 


dt^ 


dx 


^ dB^^ dB^^ dx-^ 


meaning grad^-LA^", a = 1,2, where AT = (A“^, A“^, i?“). Since grad$ and 
A^ X A^ are collinear, the vectorial expression for complete integrability condi¬ 
tions for plane non-linear PDE systems is derived: 


( 2 ) 


grad$ X (A^ x A^) = 0. 


Explicitly, this means 



= 0 


= 0 


= 0 , 


where (P, Q, i?) = A^ x A2 = {A^'^B^ - A^^B\A^^B^ - A^^B^,A^^A^^ - 
A12a 21)^ The third relation above is a consequence of the previous ones and 
the differentiation of the first and the second one with respect to B , respectively 
B is leading to the explicit integrability condition for one state plane non-linear 
PDE systems: 


(3) 


dB\ dp) dpy^ dp)' 


2.2 Multi-state variables plane PDE systems 

In the following, the main interest is to extend the previous result for linear 
PDFs having n state variables. Let Aj = (A“^ : D x P” x > R)aA=T 2 
B = (P“ : D X P" X P^ — > R)a=T 2 denote some 2x2, respectively 2x1 matrix 
fields, where i = l,n denotes the state variables index. Let us start with the 
following guasi-linear PDE system: 

(4) Ai{t,x(t),u{t))gr&dx''{t) = B{t,x{t),u{f)). 

The basic and also an original idea used here consists in separating the 
gradients of the two state variables via some additional control. More precisely. 
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it consists in rewriting the above PDEs as a split non-linear PDE system,'. 


( Ai{t,x{t),u{t))gia.dx^{t) = Vi{t), i = l,n — 1 (no sum); 
\ Anit,xit),u(t))gradx"(t) = B{t,x(t),u{t)) 


Then, by applying the result obtained in the previous Section, the complete 
integrability conditions for n-state variables quasi-linear PDE systems may be 
phrased: 


( 6 ) 


A 

dt^ 



dP 


Q 

b/* a: ^^2 


i = l,n, 


where 


{P,,Q,,R,) = {Afv^-AfvlA: 




A\^vf,det Ai), Vz = 1, n — 1 


p„ = (b^ - E“-‘ „!(()) - (s' - Ew «;(*)) i 

= Al‘ (b- - E”.-,' -Hi)) - iiH (b" - E”.',' -Hi)); 

Rji — det Aji. 


2.3 The PDE system of stress tensor in perfect plastic 
medium 


As mentioned at the very beginning of this section, the content of the paper is 
related to the stress tensor for deformations in perfect plastic plane medium. 
In this paragraph, a PDE system associated to this geometric object is defined 
and rewritten in a more manageable way and, by applying the theoretical results 
obtained above, the corresponding complete integrability conditions are phrased. 
For this, it makes total use of the inspired idea of splitting linear PDE systems 
over gradients; this separation will also prove to be of major utility in the further 
development of the paper. 

There is no novelty that the stress tensor components for perfect plastic plane 
medium are described by a constrained non-linear PDE system (equilibrium 
condition, see m, i): 


(7) 




Ido'xx d)(7xy _ 

dx dy 

dcxy dayy _ 

dx dy ~ ' 

{CTyy - axxf + daly 


4a:^ 


where Uxx, <Xyy (normal stresses) and Uxy (shear stress) denote the components of 
the symmetric stress tensor and K is the radius of the Mohr’s circle (depending, 
as Mohr’s equation proves it, on the average of the normal stresses). 

Using a classic and natural change for state variables (meaning a rotation of 
angle ^ in order to replace the given coordinate system with principal directions) 
it leads to Uxx = p — K cos (/?, Uyy = p -\- K cos (/?, Oxy = K sin p and may be 
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written the polar non-linear PDE system: 


( 8 ) 


dp OK . dip dK , n 

-r-cos ip K smip— —h sin(^ + A cos p— = 0, 

dx dx dx dy dy 

dp dK . dp dK , dp 

——h cos p — K sin 1 ^——h sin p K cos iy9-— = 0. 

dy dy dy dx dx 


It follows that {x,y) may be identified with the bi-time variable, {p,p,K) 
is the three-dimensional state variable and the matrices of the non-linear PDE 
system are 


Ai = 


1 

0 



— cos p sin p 
sin p cos p 


^3 


K sin p K cos p 
K cos p —K sin p 


It is obvious that initial control variables are missing, but some natural 
control variables emerge as consequence of the canonical split of the system 
over gradients. Applying this technique here leads to 


(9) 


dp dp 

dK dK 

-—cosp-\- —sinp = p, 

dK^ dK^ 

sm p -\- cos p = n, 

dx dy 

■ dp dp 

K smw--—h K cosp-^ = —u — u, 
dx dy 

dp dp 

K COSif— - K SlIK^-;:— = —V — ly, 

ox dy 


or, equivalent 


dp 

OX oy 


dx 

dK 


dK 


= — y^cos^J-I-i^sintp, -—= psmpi/cosp, 
dx dy 

dp 

K = —(u-\-p) sin p — cos p, 

K= —{u-\-p) COS p (v v) sin p. 
dy 


Finally, the integrability conditions for perfect plastic medium problem are the 
result of relations (6) and have the simplest expression 
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dv 


< 


— (—/iCOS(^ + sin ip) = —{psinp + vcosp); 
oy ox 


(^{u + fj,) sin p + {v + iy) cos 


dK dp 
dx dy' 


dK dp 
dy dx 



cos p — {v + v) sin p 


Remark 2.1 As direct consequence for the above integrability conditions, for 
each predefined state p = p{x,y), the rest of the state variables p and K are 
solutions for the complete integrable PDE system 


dp d ^ ■ \ 

dp d d 


3 Optimal control problems constrained 
by non-linear PDE systems 

The aim of this section is to use arbitrary non-linear PDE systems as constraints 
for optimizing cost functionals defined as multiple integrals. In the most general 
framework, the main interest consists in finding 


(11) max/(«(•)) = / X{t,x{t),u{t))dv + / g{t,x{t))da, 

“(■) Jn J^=dn 

constrained by 

_ 

Ai°‘{iix{f),u{f)) — {f) = B'^{t,x{t),u{t)), = 1,2. 

In the spirit of the breakthrough in the first Section, it is of major utility to 
replace the previous constraints with a split n-state variable quasi-linear PDE 
system 


J Ai{t,x{f),u{t))gmAx'^{f) = Vi{f), i = 1,n — 1 (no sum); 
\ A„{t,xit),u(t))gradx^(t) = B{t,x(t),u{t)) - 


The complete integrability conditions (6) define the set of admissible controls 
(with an initial component u = {u^, and some additional ones Vi = 

i = l,n-l): 


U = = {u,v) : n 


R 


N-\-2n-2 


u constrained 


by {CIC 2 )} 
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Remark 3.1 Any variational problem 

max/(a;(-)) = / X{t,x{t))dv + / g(t, x{t))da, 

^(•) Jn J-s=dn 

constrained by 

pt^i _ _ 

A^°‘{t,x{t)) — {t) = B^{t,x{t)), a,l3 = 1,2, i = l,n 

becomes, after the natural split of the PDE system and the natural addition of 
canonical control variables, an optimal control problem (with no initial controls). 
This fact was anticipated in many papers describing extensions from calculus of 
variations to optimal control (see ® m, iw, iBi, Nevertheless, 

the originality of this paper consists in the fact that explains, for the first time, 
the natural process of inserting canonical control variables compatible with the 
constraints. 


In order to solve the optimal control problem and following the classical pattern, 
some Lagrange multipliers p’‘(t) = i = l,n are introduced in 

order to define the Hamiltonian 


n—1 


H{t,x{t),u{t),p(t)) = X{t,x{t),u{t))+ ^p^p{t)v'^{t) 

+ x(t),u{t)) - ^ uf (t)^ . 

The initial variational problem is changed into an optimal control one, 
namely 

(13) max/(«(•)) = / C{t,x{t),u{t),p{t))dv + / g{t,x{t))da, 

“(■) Jn Js 


or 


max/(u(-)) = / 

«(•) Jn[ 


TT i A0ad^ 


dv+ g{t,x)da. 


Although several more general approaches on control theory have been ana¬ 
lyzed so far (see 0 , [in]), this paper bases on the assumption that the optimal 
problem ((12), (13)) admits a continuous optimal control u*(t) = {u*(t), v*{t)) G 
Int{U) which generates an optimal state x*{t). Let us consider a variation of the 
control u^{t) = u*(t) + S,h(t), where h is an arbitrary continuous vector. Since 
u*{t) G Int(U) and a continuous function over a compact set Vt is bounded, 
there exists ih> ^ such that u^{t) =G IntU{t), V|^| < This ^ is used in our 
arguments. Furthermore, let a;j(t) be the corresponding variation of the optimal 
dx^ 


state and y{t) = 
define the function: 




be the variation vector field. For |^| < let us 


HO = / At,xl{t),u^{t),p{t))dv + / g{t,x\{t))da = 
JQ Jd 


= In l^(i,xl(t),uf(t),p(t)) -p'p{t)AA{t,x\{t),u^{t))^{t) 

+git, x^{t))da. 


dv+ 
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Differentiating the function I (^) leads to 


f I dH dxi dH 


dul 

W 


d‘^ 

du'i 

+^:^{t,x^,u^) — 


du' 

For ^ = 0, it follows 
/'(O) = 




dA^°‘ dxi 

-a^rX^.ru,)^ 


r do dxi 


in 
[dA^i^ 


^^-it,x*,u*,p)y’- + ^{t,x* ,u* ,p)h°^ - ppA^°‘{t,x* ,u*)^^ 


dx 


dt° 


' dA^°‘ n dA 




0a dy" _ d 


By substituting the term —with — 
the relation above may be rewritten 




dP 


MAt)y\ 


/'(O) = + 


- vlA^^y* 


dH ^ dH 5 r , , 

+ 7)7^ 


d^y -P0^^ y 

6 

dx^ ^ ^ i9u“ ’ ” ' dt 

{p^pAt) 


n 


da 

/3a 


i i/'dAf°‘ ■ dA^^ \(Jx" I 


dx^ 


dH d 


, dA^^ dxA 

dx^ dt^^ 


dH idAf^^dx" 
du°- du°- dt°‘ 


h^dv 


+ 


/S L 


dy i A 0a 

^.-PpA. na 


y'^da. 


Since ^ = 0 is the critical point of /(C), it follows /'(O) = 0. Defining the set 
V = {p*{t)} of optimal costate variables as the set of solutions for the Cauchy 
problem 


^it,x*,u*,p*) + ^ (p»Af“ {t,x*,u*))- p>;it) ^ it,x*,u*)^it)=0, 


d 


dAl 


dx^* 


dx^ 

■^it,x*) -p^i^*it)Af°‘it,x*,u*)nait) 


dx^ ’ dt° 

= 0, Vz = l,n ( i.e. no sum over z), 


leads to the optimum critical point condition 

dn d-r*'^ 

— (t,X ,U ,p)-pp it)^^it,X ,U )— = 0 

According to the above statements, the following result may be phrased: 

Theorem 3.1 Ifu*it) is an optimal solution for the problem ((12), (13)) and 
x*{t) is the corresponding optimal state, then there exist the costate vectors 



p* = ip’j^)0-Y2 such that {x*,u*,p*) satisfies 

dx 

I Vf =Ar{t,x\u-, 

(14) 


v^*=Atit,x*,u*)^ 0 


no sum over i 


0'^ 


BP{t,x\u*) - E”=”i^r = 


dx^ 


(15) 

f)Tj f) / \ - 

— {t,x*,u*,p*) + ^{pfA^^it,x*,u*)) -p>;{t)-^{t,x\u*) 

(no sum over i), 


. dxA 
dt°‘ 


= 0 


(16) 


dP V * dB / dB I ’ 


(17) 


— {t,x*,u*,p*)-p}*{t)-^it,x*,u*)i^{t) = 0, \fa=l,N + 2n-2-, 


du°- 


,dx* 


(18) 




V / I * —*\ 

(t,X ,U jUc 


0 (no sum over i). 


Remark 3.2 For a variational problem as the one described in Remark 2.1, 
the relations above gain a much simple expression and they are much more easy 
to work with, due to the absence of initial control variables. More precisely, u 
consists only in the canonic additional control variable v. Then, 


(19) 




dx* 


(no sum over i)] 
(9a;*" 


B^{t,x*)-j::-,\t = A^„‘^{t,x*) 


dt° 


( 20 ) 


OH 




dA^^x\v*y) + —Ar{t,x*)+fg 




dx^ 


{t,x*)^— = 0 (i 
^ ' dt‘^ ( 


no sum over i 


( 21 ) 


dH 


{t,x*,v*,p*) = Q, 


Vo; = 1,2; 


( 22 ) 


dg 

(9a;® 


{t,x*) - p'p A'^F{t,x*)n, 


0 (no sum over i). 
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4 An optimal control problem constrained 
by perfect plane medium PDEs 

This section combines the results related to integrability of the PDEs describing 
the stress tensor components in perfect plane medium (Section 2) with the max¬ 
imum principle derived in Section 3, in order to completely determine a solution 
for a given optimization problem. Similar ideas have successfully been applied 
before in connection with stochastic theory, differential geometry, deformation 
theory, electromagnetic fields m, i, [13], [16], [18], [E]). It is emphasized 
once more, if necessary, the high applicability of the optimal control techniques 
in mechanical processes. 

Example 4.1 Let Q, = D 2 denote the unit disc in and let E = S '2 be its 
boundary (the unit circle). Let us consider the following variational problem: 

(23) max / (p{x,y)da, 

(V(')) JS2 


constrained by PDE system: 


(24) 


dp dK , dip dK . dp 

- 7 : -cos p + K smp— -h sm -I- A cos p-— =^, 

dx dx dx dy dy 

dp dK , dp dK . dp 

— -I- cos p — K siTvp— — h smp + K cos p-— = 0. 

dy dy dy dx dx 


The natural approach of this problem stands on the separation of the state 
gradients, via some additional, canonic control variables. Denoting these canonic 
control variables with u{x,y), v{x,y), p,{x,y) and i/{x,y), the variational prob¬ 
lem may be reconsidered as an optimal control problem constrained by: 


dp dp 

dK dK 

= —p cos p + sm p, = p sm p v cos p, 

dx dy 

dip 

K = —{u + p) sin p — (v + ly) cos p, 

K-tP = —(u + p) cos p + (v + ly) sin p. 
dy 


Multiplying each relation above with a corresponding Lagrange multiplier 
Pi{x,y), P2ix,y), ri{x,y), r2{x,y), qi{x,y), q2{x,y) respectively, the corre¬ 
sponding Hamiltonian writes 


H = piu + P 2 V + ri(^ — p cos p + ly sin p'^ + r 2 (^p sin p + v cos p'^ 

+ 9 i ( ~ + m) sin p — (v + ly) cos p'^ + q2(^ — (u + p) cos p + {v + v) sin ipj 


Application of the maximum principle for variational problems (Remark 2.2) 
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generates the following PDE systems: 


(25) 


d{Kq, 

dx 


dpi dp2 

— 0 - 

dx dy 


dri dr 2 

dip 

dx dy 

= ^^Tx 

d{Kq2) 

= sin p ^ 

dy 


' dy 


— cost^^ri^ + r 2 fJ- — qi{u + y.) + q 2 {v + v'^ , 


(26) 

and 

(27) 


Pi = qi sin ip + q2 cos (p; P2 = qi cos (p — q2 sin (p; 
ri = -q 2 ', r 2 = qi- 


dg „ 

Pini + P2n2 = ^ + P'^y = 

dg 1 V 

qini+q 2 n 2 = ^ ^ qix + q 2 y = 1, on E. 

op 


rim + C2n2 = 


dg 

dK 


Tix + r2y = 0 . 


Using relations (26), system (25) becomes 


(28) 


^ ^ = n- 

dx dy 

dq2 , dqi dp dp 

“"5—c g2-^ 
dx dy dx dy 

dqi dq2 dp dp 

"b-^ "b~ ~ -^92-^, 

dx dy dy dx 


with the particular admissible co-states pi{x,y) = y, p 2 {x,y) = —x, qi{x,y) 
r 2 {x,y) = ys\np{x,y) - xcosp{x,y) and q 2 {x,y) = -ri(a;,y) = ycosp{x,y) 
xsinp{x,y). Introducing them into the boundary constrains (27), leads to 

{y^ — x^) cos p + 2 xy sin p = 1 on 82 - 
Since (y^ — x^)^ + (2a;y)^ = (a;^ -b y^)^ = 1 to 82 , it follows 
cos p{x,y) = y'^ — x'^, sin p{x,y) = 2 xy on 82 ■ 

Extending this solution on D 2 , gives: 

cosp{x,y) = , sinv3(x,y) = , 




H- 


. 2xy 

hence p(x,y) = arctg ^5 - 5 -. 

y2 _ J.2 
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By introducing these in Remark 1.1 (combining integrability conditions with 
the initial PDE constraints, for a given state variable (p{x,y)), it follows that 
the state variable K needs to satisfy the second order PDE 


2(y^-x^) 


d^K 

dxdy 


2xy(^ 


d^K d^K\ 
dy'^ dx'^ ) 


+ 4:{y 


dK dK 
dx ^ dy J 


u 


= 0 . 


Look for particular solutions of type K{x,y) = A{x) + B{y) leads to 


K{x,y) 


aix"^ + y'^) + /3- + 7 - 


6 . 


In particular, 

1. if K{x,y) = a{x‘^ + y^), then p{x,y) = —2a{x'^ + y), therefore K = —\p', 

2. if K{x,y) = /3i, then p{x,y) = and K = p; 

3. similarly, for K{x,y) = 7 ^, p{x,y) = 7 ^ and K = p. 

4. Finally, if K{x,y) = 6, then p{x,y) = —(51n(a;^ +2/^)- 


5 Conclusions and later development 

The example analyzed above explains the utility of a multitime Pontryaguin 
maximum principle when dealing with mechanical phenomenons described by 
non-linear PDFs. The spectacular outcome is the fact that the solutions em¬ 
phasize some reasonable expectation related to the dependence of K on p (fact 
already known as basic feature in plastic deformations). 

Moreover, the approach described in this paper encourages us to apply sim¬ 
ilar techniques for even more sophisticated problems in physics (generated, for 
example by higher order PDFs) or even in differential geometry or other scien- 
tihc areas. 
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